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Abstract. We consider several new classes of viable vector field alternatives to the 
inflaton and quintessence scalar fields. Spatial vector fields are shown to be compatible 
with the cosmological anisotropy bounds if only slightly displaced from the potential 
minimum while dominant, or if driving an anisotropic expansion with nearly vanishing 
quadropole today. The Bianchi I model with a spatial field and an isotropic fluid is studied 
as a dynamical system, and several types of scaling solutions are found. On the other 
hand, time-like fields are automatically compatible with large-scale isotropy. We show 
that they can be dynamically important if non-minimal gravity couplings are taken into 
account. As an example, we reconstruct a vector-Gauss-Bonnet model which generates 
the concordance model acceleration at late times and supports an inflationary epoch at 
high curvatures. The evolution of the vortical perturbations in these models is computed. 
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1. Introduction 

Most inflationary and dark energy models [1] are based on scalar fields [2]. They, as 
conventional cosmological sources, are fully characterized by their equation of state and 
speed of sound. For a more general phenomenology, one should also take into account 
possible imperfect properties of a fluid [3, 4]. Indeed, it has been proposed that dark energy 
might have such a characteristic [5] , and the possibilities to distinguish it observationally 
have been examined recently [6,7]. In Ref.[8] we contemplated on the possibility of such 
a significant anisotropy: that the expansion rate would be better described as a direction- 
dependent quantity within the Bianchi I model. The phenomenological set-up there 
could describe viscous fluids, Yang-Mills fields, strings or cosmology with noncommutative 
infrared properties. It was found that significant anisotropy may be allowed in some cases, 
and that these models have potential to explain the anomalous cosmological observations 
[9,10]. In the future, one could optimistically expect to detect signatures of possible 
early anisotropies in the polarization, and of possible late anisotropies in the supernova 
luminosity distance-relationship data. 

In the present study we take a complementary approach and rather than working 
within a framework of phenomenological parametrization, we study vector field theories 
with explicitly written Lagrangians as our starting point. In particular, we look for 1) 
solutions which may be compatible with symmetries of cosmology, and then 2) which of 
these solutions might be used to model infiation or dark energy. The resulting models 
can then be divided into two classes, those with space-hke and those with time-like field 
components. Both cases then seem to violate Lorentz invariance. However, this violation 
is already done when choosing the cosmological background, not when adding the vector. 
The cosmological metric picks up a direction for time. The vector fields we consider are 
then particular solutions of a Lorentz-invariant theory compatible with the broken space- 
time symmetry of the metric. The conceptual difference with the Aether [11] and such 
theories arises from the unit- norm requirement, which is achieved by adding a norm- fixing 
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term into the Lagrangian that then forces the field to be always time-like, even in vacuum. 
In the present study we do not impose the unit-norm requirement. See [12-14] about 
spontaneous Lorentz violation in string theory and gravity. 

Space-hke vector fields in cosmology have been used to model infiation by Ford and 
others [15,16]. We will therefore focus on the present study more on the dark energy 
era, where the qualitative differences are that matter cannot be neglected and that the 
cosmologies do not have to isotropize. For completeness we also include a coupling to 
matter. In particular, we work within the Bianchi I model, and look for scaling solutions. 
This approach differs from previous studies which have been confined to the case of the 
exact Friedmann-Lemaitre-Robertson- Walker (FLRW) symmetry. Space-like fields are 
compatible with this symmetry only in the case of a so called triad, which has three 
identical vector fields for each spatial direction [17]. This can be used to model dark 
energy [17-19]. Inflation has been also considered in the case that a large amount of 
randomly oriented vector flelds result in a average FLRW spacetime [20]. We note also 
that a vector fleld with quintessence can be considered [21], and that adding nonlinear 
terms to the electrodynamic Lagrangian could cause the universe to accelerate [22]. 

Time-like fields have acquired a lot of interest in cosmology [23-31]. If the field is 
canonical and minimally coupled, the equation of motion for its time-component becomes 
trivial. Therefore one has to cither couple the field non-minimally or add non-canonical 
terms (these modifications are equivalent in some cases). Within modified gravity models 
attempting to eliminate dark matter, vector field components have been promoted to a 
central role due to their ability to mimic the observed gravitational lensing phenomena 
from dark matter [25, 32, 33]. 

Recently non-minimally coupled time-like vector fields have been taken under 
consideration as dark energy candidates [34-36]. Even if the vector field does not drive 
the infiationary expansion, it might have a role in slowing the expansion rate or on 
the generating perturbations [23,37-40]. Recently the observed anomalies in the CMB 
have inspired a number of studies about perturbations in an anisotropic background [41- 
44], generation of structure during anisotropic infiation [45-51] and on the formation of 
structure in the presence of anisotropic sources [52,53]. In the present study we focus on 
role of the vector fields as driving the background expansion, especially the accelerating 
expansion associated with dark energy and inflation. The evolution of the spin-1 type 
perturbation in such a case is determined and we find these modes are decaying. 

We study cases where the field is coupled to gravity non-minimally. These vector- 
tensor couplings can cause the field to effectively break the usual energy conditions. 
Therefore they could be applied to model non-singular cosmologies and to construct worm- 
holes [54]. If the photon has such couplings during infiation, large-scale magnetic fields 
could be generated [55-57]. In the present universe, the vector-tensor theories can be 
severely constrained by Solar system experiments [58,59]. We consider in particular two 
classes of these theories: with a couphng to the Ricci scalar and to the Gauss-Bonnet 
invariant. The latter has special interest, since the Gauss-Bonnet term has theoretical 
and phenomenologically desirable properties. In particular, its presence in the action does 
not immediately imply existence of ghosts and severe conflicts with observational data on 



CONTENTS 



4 



gravitation. Our results include that this model could generate both the inflation and the 
present acceleration of the universe. The outline of the paper is simple: in section 2 we 
derive our results concerning space-like models, in section 3 those about time-like models, 
and in section 4 we conclude. 



2. Space-like Fields 

2.1. Minimally Coupled Vectors 
Consider the vector field action 



(1) 



IGTrG 4 

where F^j, — d^A^, — di^A/j, and A"^ — A^A^^. The field equations are 

G,. = SttG (T- + (2) 

where the energy-momentum tensor of the vector field follows by varying with respect to 
the metric 

T^, = F^F;; + 2V\A^)A,A, - + V{A^)^ g,.. (3) 

The equations of motion for the four vector components can be written as 

V^F'^'^ - 2A''V\A^) = 0. (4) 

These follow either by setting the divergence of (3) to zero or setting the variation of (1) 
with respect to the field to zero. 

We want to study the consequences of such matter sources in a Bianchi I type universe. 
Such is a generahzation of the FLRW case with three different expansion rates. Thus we 
allow for a possible anisotropy. Then the metric may be written as 

ds'' = -de + a''{t)dx'' + + ^{t)dz' (5) 

3 

= - dt" + Y^{aiXif , (6) 

i=l 

where we have introduced an obvious labeling of the spatial coordinates and the metric 
components with Latin letters. Following this convention, we denote A^ = {(f), A^). Then 
the components of the energy-momentum tensor may then be written as 

^o'o = I E i A' + ^(^') + ^V'(A")<I>"^ (7) 

1=1 ' 

T,i = 2V'{A')cl>A,, (8) 

T,j= - + 2V'{A')AAj + a,a, ' ^(^V (9) 

\^ k=i "fe / 

Because the symmetries of the metric (6) do not allow a velocity field Toj, one immediately 

sees that we are restricted to consider either a purely time-like vector, A^ = (0,0), or 

a space-like vector with = 0. However, the former case reduces to triviality since the 
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equation of motion (4) for the nonzero vector component states that V'{A^)(j) = (since 
one has V^F'^*^ = 0). Therefore we will consider the space-like vector fields. For them, 
the off-diagonal components of Tik should vanish. In principle, one could then consider 
nontrivial (space-like) vector fields which satisfy 

-AiAj + 2V'{A^)AiAj = 0. 

Then the evolution of the vector field, Ai{t), would be determined a priori without reference 
to the other matter (which then might or might not be consistently added into the system) . 
Such possibilities seem contrived and we exclude them from our present considerations. 
Thus, we have found that only space-like vector fields parallel to one of the principal axis 
could be allowed. For the FLRW metric all the diagonal components of T^j should be 
equal, which then requires three vector fields, one in each coordinate direction, of exactly 
equal magnitudes, as the triads of Armendariz-Picon [17]. With the metric Eq.(6), the 
only constraint is that every vector field should be along a coordinate axis. Any system 
of these vector fields is thus required to have an anisotropic equation of state, unless it 
reduces to the very special triad case. Let us consider the case 741 = ^43 = 0, A2 — B. We 
then define the additional dimensionless parameters 

where H is the average expansion rate 

It is also useful to define the fractional difference between expansion rates as 

b 




b 



(12) 



This is a simple quantification of anisotropy with immediate generalization to the invariant 
shear. We can then make contact with the parametrization of [8, 60] by noting that 

X'-2Y -X' + Y, 

where w is the equation of state along x-direction, and 6 is the difference of the equations of 
state along y and x directions. Generalization to multiple vector fields is straightforward. 



2.2. Scaling Solutions with and without Matter Couplings 

In this subsection we consider space-like cosmological vector fields. Similar anisotropic 
infiation has been considered in the early universe [15,16]. Therefore we focus on the 
present study more on the dark energy era, where the qualitative differences are that 
matter cannot be neglected and that the cosmologies do not have to isotropize. Our main 
aim is to find scahng solutions. There are several reasons to be interested scahng solutions. 
1) Simply enough, if they are attractors, they may describe the realistic dynamics of the 
field. In any case they tell us about the phase structure of a given model. 2) Scaling 
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solutions have also the possibility to alleviate the coincidence problem. If dark energy 
and dark matter had similar (even within an order of magnitude or two) energy densities 
in the past, it seems less discomforting that the energy densities are about equal just 
today. However, it may still remain unexplained why the transition to acceleration begun 
just at small redshifts - this is even if the dark matter and dark energy would continue 
scaling in the future, i.e. they would always have (roughly) equal magnitudes. Still, 
the coincidence problem a is good motivation to study scaling solutions [61-65]. 3) In 
addition, the possibility of early dark energy, which the scaling solutions automatically 
incorporate, is phenomenologically interesting [66]. A model predicting the presence of a 
dark energy component in the earlier universe can already be strongly constrained, and it 
may feature new signatures which might be useful in distinguishing between alternative 
models of acceleration. Scaling solutions are known to exist in FLRW universes for specific 
forms of scalar field Lagrangians and these have been applied in attempts to address the 
coincidence problem. Such a property could be shared by space-like vector fields in the 
Bianchi I background. For completeness we also include a coupling to matter. However, 
the presence of a coupling is not a necessary condition for this, as will be seen below. 
Consider the action 



/ 



(14) 



167rG' 4 

where we have included a general matter Lagrangian L.^ and the coupling function q{A^). 
To proceed, we write down the complete evolution equations using the variables (10). The 
amount of matter can be now described by a generalization of the usual fractional energy 
density, 

(15) 

The Priedmann equation may then be written in the form 

3 = 3fi^ + ^X2 + F + li?2. (16) 

Our time variable is chosen to be a generalization of the e-fold number. The prime will 
denote a derivative with respect to that, 

^ dlogahc ^ ^ 

It is convenient to derive an average effective equation of state, which then describes the 
overall expansion rate of the universe. 

^e// = -^§-l, (18) 

Finally, we define two dimensionless variables, 

/?^^, g ^ ^g'(/l2)/3, (19) 
which are the comoving field /3 and the coupling function q. The k, — l/(87rG). 
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The system of equations then becomes 



X' = 
Y' = 
R' = 
9!^ = 



X 



2q{l- 3Wm) ^r, 



Yi 

2 — y + — 

H 13 



X-(3{l--R 



H' 



R + X^-Yi, 



fir 



H' 2 
2— + 3(l + «;™) - -Xq\ 



2,Wm) 

Y' = 



(20) 
(21) 
(22) 
(23) 

0. It is easy to see 



For a conventional scahng solution one should have X' 
that this makes sense only if also Rl = 0. Then the derivative of the average Hubble rate 
is a constant, 2H'/H — —3(1 + Weff), which is related to the effective (average) equation 
of state for the universe, Wg//- Since also w and w + 6 now both correspond to the dark 
energy equations of state, 6 and thus Yi must be a constant as well. One deduces that this 
could be satisfied only if either 



2\n 



= 0, 2)Y ~ (A^) 



or 



3)y = Yi = 0. 



Below in table 1 we list the results for each case, but before we also derive and discuss 
these in more detail. 

• In the case 1) that the comoving field stays constant, the square bracket term in 
Eq.(21) vanishes (since, as one can readily check, the derivative of the comoving field 
is proportional to the square bracket term). One uses this fact to relate the comoving 
field to X and R. Moreover, now either Y = or H' = because Eq.(21) must be 
satisfied. The former possibility of vanishing potential goes under the case 3) which 
we consider below. The latter possibility would correspond to the case Weff — — 1, 
which in this regard resembles a de Sitter-like space. However, since Eq.(22) sets now 



R^ 



X^-Y^ 



the solution is anisotropic unless = Yi. From Eq.(23) one sees that for a scaling 
solution with Q^n 7^ 0, one would have to force a coupling to keep the average Hubble 
rate constant. Since g is a function of the comoving field, the coupling term is now 
a constant. Plugging this into Eq.(20) only gives us an equivalent of the Hubble 
equation as a consistency check. Therefore we have a whole set of solutions, which 
can be parametrized by the constants Yi and q. We summarize these solutions as 



X 
Y 



-1, 

^ 9(1 + w^) 

2^(1 - 3Wm) ' 

= |648g^ {l-3wr, 



'1 + w^' 



- 6561 {-1 + Wm) (1 + m 

(9 - 27Wm + 4:WmY,) 



mJ ~ 
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R 



- 16g^ (1 - Swrn)' {Y,' - 81) 1 } .^^.,.\. .4 . (24) 

J J 432/? (g - 3qwm) 

4q\l - 3w)2 



3 

(3(l + wJ (4g2 (1 _ swj' (Y, - 9) - 81 (1 + wmf) + (Sw^ - if Y, 



8/3 (g - Sqwmf 

The expression for the most general case is comphcated. However, one notices that 
this scahng sohition is unhkely to be able to describe a realistic era between the 
inflation and dark energy domination, since it implies that H is equal to a constant, 
and requires a coupling to matter, which can easily be problematic 
• In the power-law potential case 2) one notes that V = is consistent with Eq.(21) 
only if X = (unless we would be back into cither of the two other cases). Then 
w = —1 and 6 = 2n/3. The coupling then does not affect the matter scaling since the 
field is constant, and Eq.(23) tells that Wgff = Wm, since we would like to consider 
scaling solution with nonzero Q^. The constancy of R requires, from Eq.(22) that 

AnY 

R 



Wm-l 

Plugging this back into Eq.(21) we get that 

2n-3{l + Wm) 



Y = 9{wm - 1) 



16n2 



This should of course be positive, imposing a constraint n < 3/2 if Wm — 0. The 
amount of matter one then finds from the Priedmann Eq.(16), is given by 



8 V n 



This may be positive for Wm = when n < —(3 + -\/33)/4 ~ —2.19 or n > 
(—3 + v^33)/4 ~ 0.69. Thus we have found that there exists a scaling solution for 
a vector coupled to matter, which may be described as 



X =0, 



^ ^ 9{wm-l)[2n-3{l + wm)] . . 



R 

0*7 



-{9-6n + 9wm) 
4^1 ' 
3 {2n^ - 3(1 + Wm) + n{3 + Wm)) 



8n2 

However, since we should have X' = X = 0, one may be concerned with Eq.(20). 
We mean that the kinetic contribution to the total energy should stay negligible, and 
thus terms in the RHS Eq.(20) should somehow cancel. One option would to be to 

I A possible solution to this issue would be to use chameleon vector fields [67-69]. 
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introduce a coupling which cancels the effect of the potential that would otherwise 
set the kinetic contribution X varying: Pq = —nY/{{l — 3wm)^m)- This could be 
consistent only for a coupling of the specific form q{A'^) ~ loglA^). However, from a 
practical point of view we rather notice that Yi/j3 can be approximately neglected if 
the field has run to large enough values. We have found that this approximate scaling 
behaviour indeed arises naturally and is an attractor for the system for n < —2.2. 
These considerations can be confirmed by numerical computation. In Fig.(l) we show 
some models featuring the scaling solution Eq.(25). 
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Figure 1. Convergence to the scaling attractor Eq.(25) and the transition to the 
accelerating era for minimally coupled space-like vector fields. The potential is of the 
form V{x) = F+a;" -I- F-a:"". In the LHS figure, n = 3 and the RHS figure, n = 10. The 
solid (black) lines are the ratios U of the vector to the total density, the dashed (red) 
lines are the average total equations of state Wef / , and the dash-dotted (blue) lines are 
the expansion normalized shears R. The properties of the scaling attractor depend only 
on n. The transition to acceleration depends on the scale V+ which we have tuned such 
that the vector dominance is taking place near the present time {x ~ 0). 



• Let US finally consider the case 3) that the potential can be neglected. We see directly 
that for such a solution w = 1 and 6 = —2/3. Again Eq.(23) gives the effective 
equation of state, Weff = Wm — 2qX{l — 3wm)/Q which now may be influenced by the 
coupling. Then Eq.(22) gives 

9 + 2qX -3Wm{3 + 2qX)' 

Inserting fi^ = 1 - X'^/Q - R^/9 in Eq.(20) would now give an expression for X, but 
this is generally quite complicated. Let us first look at the case that coupling vanishes. 
Then 

Weff = Wm, 

3(3 - Wm) , . 

^im = , (26) 
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Solution 


Equations 


^eff 


w 


6 


Existence 


R 


Comoving field 


(24) 


-1 


w{P,...) 


w{f3,...) 


Needs coupling 




Potential driven 


(25) 


Wfji 


-1 


1/3 


Large field 


3 9 
2 4n 


Kinetically driven 


(26) 




1 


-2/3 


Not for dust 





Table 1. Scaling solutions for space-like vector fields and matter. The last column on the 
right shows the shear R in the dust filled universe. All of these solution exhibit anisotropy. 



3(-l + 3Wm) 
^ = 4 ' 

X = ± ^ + (4 - SWm)Wm. 

The case of immediate concern to us is a background of dust, Wm = 0, but there this 
solution does not exist since X should of course be real. Thus we are led to consider 
a case that the field is nonminimally coupled to dust, Wm = 0, q 0. This requires 
that the coupling function q is constant, which is possible, if /3 is not a constant, only 
for the special form of the coupling q{A'^) ~ ^/{A^. We again dismiss this possibihty. 

We focus on the most relevant case, where dust dominates and there are no nonminimal 
couphngs. We show that the scahng solution (25) is an attractor. For that purpose, we 
ehminate the H'/ H and Qjn terms from the system (20-23) and get 

X' ^ — (-6 -8R + 2i?2 + X^-6Y + 4nY) + 0{h, (27) 

(28) 

y = K (18 + 2R^ + X^-6Y + 4n(-3 + 2R + ¥)) + O(^), (29) 
6 p 

(30) 

R! + ^ {2R^ + 12(^2 - 2nY) + R{-18 + X^ - 6Y + 4nF)) + C(^). (31) 

We consider then areas where /3 is large and terms proportional to its inverse become 
negligible. Then consider the expansion 

X = Xex + 0{€^) = {R, X, Y) + e(x, y, z) + C>(e'), (32) 

where background marked with bars is given by the solution (25). Then X = at the 
background order, and to first order in e we get 

X = Mx, (33) 

where the three eigenvalues of the matrix M are given by 

[Ml, M2, M3] = - 1) ' - i (1 + |^V(4n-3)(n2 + n-3)^ , (34) 
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The real parts of the eigenvalues are always positive for n < — 1(3 + \/33), so the point is a 
local sink. This confirms that the scaling solution is an attractor given that the exponent 
n is negative enough. 

We have thus systematically considered every possible way towards a scaling solution 
of a cosmological (space-hke) vector field which might be coupled to matter. Several 
possibilities were noticed, and of particular relevance is the scaling solution seems to be 
Eq.(25), which is valid in the limit of large enough field values. We should however note 
that several kinds of different behaviours, anisotropic or not, with matter domination 
fim, = 1, or with vanishing matter contribution VL^ = 0, were left this time without explicit 
treatment, since we were only looking for a scaling universe. In addition, more interesting 
solutions probably exist, in which the vector field tracks the background matter (in analogy 
to the scalar field tracker quintessence) without the tracking being exact (which is precisely 
what we mean by scaling solution). We summarize the scaling solutions in table 1. 

2.3. Examples of models satisfying the anisotropy constraints 

In previous study [52] it was shown that the magnitude of the quadropole resulting from 
anisotropic background expansion is given by 

where el is the eccentricity, 

e2^(^)^l. (37) 

which tells how ellipsoidal the last scattering surface (denoted by star) appears to us (the 
scale factors normalized to unity today) . The constraint on the magnitude of quadropole is 
then Q2 ^ 2.72 10~^ However, it is in principle possible for arbitrarily anisotropic expansion 
to escape detection from CMB as long as the expansion rates evolve in such a way that e* 
remains zero. In other words, the quadrupole vanishes, if each scale factor has expanded - 
no matter how anisotropically - the same amount since the last scattering. 

This is demonstrated explicitly with examples in Fig. (2). In the left hand side 
(LHS) there is an example of an exponential potential featuring large anisotropics at late 
times after the acceleration has begun. Although the anisotropy in such case has been 
perturbatively small until today, such models could not be straightforwardly described 
within the standard approach of perturbing a FLRW metric. 

Another example, in the right hand side (RHS) Fig. (2), shows that the skewness of the 
universe, as defined by Eq.(12), could be large, while the eccentricity, as defined by Eq.(37), 
might well vanish just at the present. The latter is related to the integral log(a/6) = J Rdx, 
which might even be zero on average, if the shear R happens to oscillate suitably. This 
could naturally happen, for instance, if dark energy is a vector field slightly displaced from 
the minimum of its potential. The model depicted in the RHS of Fig. (2) exhibits this 
property and is constructed as a simple, minimally coupled and one-component vector 
field introduced in the previous section 2.1. This is a fine-tuned model, but proves that 
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there exist expansion histories exhibiting very large anisotropies with dynamics that leave 
the quadrupole amplitude, Eq.(36), nearly or completely unchanged. 




Figure 2. Two types of scenarios witli a minimally coupled vector field that satisfy the 
CMBR quadrupole constraint. The solid (black) lines are the vector field density fraction 
U ~ PA/iPm + Pa), the dashed (red) lines are the effective total equations of state of 
the universe Wef / , and the dash-dotted (blue) lines describe the evolution of eccentricity, 
E = 500e^. In the LHS figure, the displacement of the field from the minimum of its 
potential is important only at late times and only then significant anisotropies begin to 
form. The potential here is an exponential, V{x) ^ e^^. In the RHS figure, the potential 
is a double-power law V{x) ~ cc"'* -I- V+x. The dynamics of the field is such that though 
there are significant anisotropies, the eccentricity at the present is nearly zero. Note 
that we have chosen both the models in such a way that they feature a recently begun 
acceleration. Thus these models are viable though finely tuned. 

The quadrupole constraint is applied here to constrain the overall amount of 
anisotropy, and we remark that the detailed anisotropy patterns which ensue within these 
constrains remain uncovered. The general evolution equations for the overdensities in the 
presence of anisotropic stresses have been presented [52], and one may observe that when 
the anisotropy is very small (as is forced by observations) the crucial modification comes 
from the anisotropic stress of the matter source. It is this quantity which sets all the 
perturbation modes direction-dependent. This will have impact on the CMB spectra, in 
such a way that the correlations between adjacent i modes become nonzero. Also unlike 
in the usual statistically isotropic case, different m-modes do not decouple. The detailed 
computation of anisotropies from the direction-dependent inhomogeneities [52] is left for 
forthcoming studies. 

To summarize, we constructed two classes of viable space-like vector dark energy 
models, described by the features: 1) The field evolves only very recently, and thus 
the impact of anisotropy on observations is strongly suppressed, or 2) There is earlier 
anisotropy, but in such a way that it's effect cancels today. Both cases require some fine 
tuning. Note however, that the quadrupole constraint suffices to estimate the amount of 
fine tuning. The reason is that only the quadropole is of the order of the background 
anisotropy in Bianchi I models. If the quadropole is not too large, the higher multipoles 
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are of the observed order 10~^ [52]. However, at the few lowest multipoles at least, there 
can appear some statistically anisotropic features. Thus we have found explicit viable 
mechanisms to generate such anomalies in the late universe. 

For completeness, one may add that an approximate isotropy may also be achieved by 
3) considering "a triad" of identical fields [17] or possibly 4) a huge number of randomly 
oriented fields [20]. In addition, let us note the case 1) was mentioned already by Ford 
in the context of inflation: a very slowly rolling field resembles of course a cosmological 
constant. 



3. Time-like Fields 



3.1. Vector Field Coupled to the Curvature 

One might also allow nonminimal couplings of the vector and gravity. Couplings of a 
massless Maxwell vector field have been recently considered in an f{R) model that could 
feature both the inflation and late-time acceleration of the universe, where the non-minimal 
couphng of the Maxwell field could then generate large-scale magnetic fields [56, 70]. Here 
we allow also a potential for the vector field, and do not couple the kinetic term of the field 
but add an interaction with the Ricci scalar it! and the field A., 



d^x^J—g 



1/1 , .oA „ 1 



+ uj{A') R - -F^^Ff^" - V{A') + L 



(38) 



2 VSttG ' 7 4 

This is a special case of a vector-tensor theory [11,25] satisfying the following three 
conditions§ 1) the Lagrangian density is a four-scalar 2) the resulting theory is metric 3) 
there are no higher than second derivatives in the resulting field equations [58]. The most 
general theory would allow also a coupling of the form A^^A^R^v\\. In the next subsection 
we however present a new model, which is of a quadratic form in the curvature invariants, 
but also seems to satisfy these requirements (depending on how to interpret condition 2)). 

The contribution of the coupling term uj to the field equations can be presented as an 
effective energy-momentum tensor, 

G,. = 8nG (T- + T^t + T;.) , (39) 

where is given by Eq.(7) and T^^ reads 

= -uG^, + (V^V, - g^^a) u - JA^A^R. (40) 

To write the components explicitly for the line element (6), it is useful to introduce the 
notation ^4^ = (0, OjAj) for the field. We then find 

To^o= - 6iJ (^A ■ A - a;' - ^ (9i/^ - H ■ H) - Vi?, (41) 
= - (t>AiJR, (42) 

§ Whether the additional condition that the free-field energies are positive for both the metric and the 
vector can imposes further constraints on the form of uj. 

II Which in fact equals the kinetic term A^^A^Rfj^^ (Va^")^ — Vc^'^V/jA", as is seen by using the 
geometric identity [V^, Vj.]A" = R^/^ij^^A^ after a partial integration of the action. 
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T^j = 2aiaj (^A • A - 0^) (3ii" - i/^) + A • A + A • A - 0^ - 00 u'Sij (43) 
- aittj (3H + ^H'^ + hl-H- Hi- SHH^ 6iju (44) 

+ 4a^ (a ■ A^ 5iju" - AiAjCu'R, 

where R — 9if^ + 6H + H • H, and we have used the obvious notation for the 3- vectors H, 
H = {Hi, H2, H3) and similarly for A. 

The equation of motion for the time component of the field is 

(l){2V'{A'')-u;{A'')R) =0. (45) 

The Goi component of the Einstein equations is identically satisfied if either the space-like 
or the time-like vector components vanish. The condition Gij = then gives 

-AiAj + {2V\A'^) - uj{A'^)R) AiAj = 0. (46) 

which translates into condition AiAj = 0, which should be satisfied for all i ^ j- Given 
this condition and = 0, an arbitrary vector field is compatible with the Bianchi I metric. 
These requirements are identically satisfied if we have only one evolving spatial vector 
component. This excludes the FLRW metric. 

In the remainder of this section we consider the special case of FLRW universe 
including a (solely) time-like field. Similar vector models have been recently considered by 
Bohmer and Harko [34] to successfully model the late acceleration of the universe while 
satisfying the the stringent Solar system constraints. Recently Jimenez and Maroto [36] 
have also noted that a vector-like dark energy could avoid some of the fine-tuning problems 
present in scalar field quintessence type models. They showed a time-like vector field with 
a non-standard kinetic term could explain the SNIa data even without a potential. 

A time-like field does break the Lorentz invariance as well as the spatial fields 
considered in the previous subsection, since picking up the time coordinate introduces 
a preferred frame. However, in an FLRW universe one can treat the vector as an isotropic 
field whose time component resembles a nonstandard scalar field (there will be interesting 
changes when perturbations are taken into account). 

Summing the contributions from Eqs.(7) and (41) the energy density is 

= 6//00 - SH'^u + V, (47) 

where we have used the equation of motion 2V' — cu'R to eliminate the last term in 
Eq.(41). One notes that this solution does not exist for the case of a massive vector field 
V{x) = ^m^x linearly coupled to the curvature, ui{x) = ujox except in the special case of 
constant curvature. However, the more general case of nonlinear coupling and/or nonlinear 
potential these kind of dynamical solutions do exist. The pressure of the field follows from 
Eqs.(9) and (43), 

p^^-V-2 (2H^(j) - 0^ - uj' + 4: (^00) \"+ (2^ + H^^ uj. (48) 
One may check that the sum + satisfies the consistency relation 
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which follows from the covariant conservation of matter that applies now since the is 
minimally coupled [71]. An interpretation is that the vector only changes the geometry. 
The conservation law may be written as 

p^ + p^^u + H^i^-^y . (49) 

Let us first first look at de Sitter solutions assuming we may neglect p,n. Wc then set 
the Hubble parameter to a constant, Hq. Firstly, one notes that an equation of motion is 
identically satisfied if the field is constant. Then the coupling follows from (49) inserted 
into the Friedmann equation as 

= Vl{?>Hl) - (87rG)-\ 

Then, consistently, Eq.(48) gives = —SHq/ISttG). This holds independently of the 
potential. One is thus able to mimic a cosmological constant with a time-like vector 
coupling, without having to have any potential present. Secondly, one may also mimic a 
constant term even if the field is rolling. Then, however, the equation of motion gives the 
nontrivial constraint V — GHqCu' and one needs a potential 

V = GH^uj + Vq. 

If we for instance assume that the coupling has a power-law form, 

a;(0) = llIi4>'^ + uJo, 

then Eq.(49) tells us the field grows (or decreases, if o; ^ —1) exponentially with time, 
(p ~ g-Hot/a^ From Eq.(47) one may then deduce that if uji — a/2 and that if uq — 0, then 
Vq 7^ 0. Similar results have been obtained in the unit-norm scalar-tensor-vector inflation 
model studied in Ref.[38], where it was noticed that a vector energy may stay constant 
even if the field is rolling down a potential, and on the other hand that such field can drive 
inflation even without a potential. 

Finally we will turn to late time universe where pm should be taken into account too. 
We look for a model which reproduces the background expansion of the ACDM model. 
The sum (49) should again vanish. The derivative of the Hubble rate would be given by 
H — |(A — 3i?^), where A = A/(87rG'). Considering then the coupling as a function of the 
Hubble rate, cu — oj{H), we find that its evolution equation: 

^(A - 3H^)uj"{H) - AHCj'{H) + 2uj{H) = (50) 

is of the Legendre form and thus has solutions in terms of the Legendre functions. 
Therefore, there seems to be no simple form for the curvature-coupling function of time-like 
vector field mimicking the ACDM expansion. From the (Legendre type) solution of Eq.(50) 
one could then deduce, using the equation of motion, 2V' — cu'R, the form of the potential. 
But since these analytic forms would be unappealing, we don't write them down here. 
The bottom line is that there exist solutions with nontrivial vector field dynamics, with 
the energy density Eq.(47) residing in the field which however remains a constant. The 
model can therefore share an identical background with the ACDM model at late times in 
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addition to possibly driving an inflationary period at early times. In the following we will 
see, that with another kind of a curvature coupling, this can be achieved with polynomial 
forms for the coupling and potential functions. Such might appear more naturally in the 
low energy corrections to the gravitational action. We leave the possible generation of 
magnetic fields in this model into future studies. 



3.2. A Vector - Gauss Bonnet Model 



The Gauss- Bonnet topological invariant Rq^ = R^'^^'^ R^i^pcr — 4:R^''Ri^i, + R^ can have 
dynamical consequences in four dimensions if it interacts with the matter sector or has 
self-interactions. Dark energy cosmologies with self-interactions of F^q^ and with scalar 
couplings have been considered previously [26, 72-74] . Here we introduce a possible 
coupling to a vector and write the action as 



-9 



1 



R - f{A')R',,s -F.^F^^'' - V{A') + 



(51) 



In analogy with the discussions in the previous subsections, we now find an effective energy 
momentum tensor 



TL = - 8 



1 



R 



(52) 



- 4 [2/;„^i?"^ - UfR] g^, + 2fA^A,RlB. 

To case notation, we define g* = —abc. The determinant of the metric is in our Bianchi 
I background g = —abc. We then have Rq^ = Sg*/g. The components of the tensor (52) 



can be written as 



rnf 
-'-00 



TL 



48/'— 
9 

Wf^Aiai 
9 



A- A 



b) + 16/'^^ 
/ a 



(53) 
(54) 



T/- ^ - IQaiaj 



2 9 



(A.A-0^)-y + 2r|^ (a.A-#) 



5., 



9 



(55) 



Again the coexistence of space and time vector components is not generally possible, but 
in Bianchi I spacetimes one may allow one space component to evolve simultaneously with 
the vector field. 

Here we specialize to what may be the simplest possible model, a time-like vector field 
in FLRW spacetime. Picking up the relevant terms from Eq.(7) and Eq.(53), we get an 
energy density 

= -48i/3#/' + V, (56) 

where we have used the equation of motion V — —f'R^B drop a couple of terms. 
Again one might note that this solution does not exist for the case of a massive vector 
field V{x) = ^m^x linearly coupled to the Gauss-Bonnet invariant, f{x) = fgx except in 
the special case of constant Rqb- However, the more general case of nonlinear coupling 



CONTENTS 



17 



and/or nonlinear potential these kind of dynamical solutions do exist. The pressure of the 
field can be read from Eqs. (9) and (55), 

= 16 [h + H^'^ 00 + H\4)'' + #)] /' - 32H\^(j>ff" - V. (57) 
One can check that the sum 

p^+p^ = 16H (2H-H^^f-^Hf (58) 
satisfies the consistency relation 

P4> + P<P = -P4'/(^H). 

In the remainder of this subsection wc will attempt reconstructions of a coupling 
and the potential in our action (51). For reconstruction schemes in other string-inspired 
cosmologies taking into account (scalar) couplings of Rq^ term see Refs. [75-77]. de 
Sitter solutions in a vacuum would be relevant for the early inflation. Now it is clear that 
if the field is constant, one would need a potential to achieve infiation. This is because 
as a topological invariant, the Gauss- Bonnet term cannot contribute to the dynamics if 
the couphng is a constant. However, if the field rolls down its potential, the coupling 
affects the dynamics, though possibly having a constant energy density with potential and 
thus yielding a de Sitter space as a solution, analogously to the Ricci-coupled case of the 
previous subsection. Let us pursue this first. 

The vanishing of the sum (58) requires 

m = /o + /ie^^°*, 

where /o, /i and Hq are constants. Since /o is classically irrelevant, we drop it. Then 
integrating the equation of motion gives us 



V{t) ^Vo- 2AH^he 



-2Hot 



Here the constant Vq corresponds effectively to a cosmological constant and we drop it 
also. If then we insert our result into the Priedmann equation, Hq — SirGp^/S, where the 
energy density is given by Eq.(56), we obtain an equation for the field as a function of 
time. The result is 

We then find that the potential and the coupling may have a simple inverse form, 

^^^'^"Lg(1-4|A|)' ^^^'^^ 6i7TGHi{4\A\-iy 
There thus exists exact de Sitter solutions when the field is evolving. In fact the field is 
exponentially increasing as a function of the cosmic time t. Therefore additional terms, 
say positive powers of \A\, could naturally terminate the de Sitter phase and end infiation. 
Hence it could originate a possible graceful exit from infiation. 

Next we aim to find a vector-Gauss Bonnet model which produces, at late times, 
exactly the same background evolution as the ACDM model. Since the vector then mimics 



CONTENTS 



18 



a constant, the sum (58) should vanish. We want to have H = [C — 3H^), where the 
constant C = 8nX is proportional to the corresponding cosmological constant A. If we 
again think of the coupling as a function of the Hubble rate / = g{H), we can solve it as 

{C-AH')g{H)-\{C-3H')g\H)^0 



^=4 

g V C- 3H^ 

2\2/3 



dH 
If 



^ g{H) ^H\C- 3H^) 
Furthermore, using the ACDM Hubble rate 



(61) 



3H = {3H^ - C)e-^^ + C, 



(C - 3H^y'^' - (l + Ce-^^) 



integrating once more we get the coupling as a function of x: 

=^ f(x) = /o + /i (220e-2^ + 264(76-5^ + 165(7^6-^^ + d^e"^^^) , (62) 

where /o and /i are integration constants and C = 3Hq/C — 1. The constant /o is 
dynamically irrelevant and we drop it. Our model has two free functions: the potential V 
and the couphng /. The solution Eq.(62) tells how the couphng has to evolve as the scale 
factor expands like in ACDM cosmology. How this corresponds to the evolution of the field 
0, depends on the functional form of f{A'^)- We will choose its form as a polynomial in 
such a way that (f) ~ e~^/^: 

f{A'') ~ [220(-^') + 264(7(-A2)5/2 + IQdC^-A'y + C^-A^^/^j . (63) 

Then, from the equation of motion V = f'G we find that 

V\A'') ~ [-2 + C{-A''f/^] \l + C{-Ay/^] [40 + 120C'(-A2)3/2 

+ 120C\-A^f + C^{-AY^^j , (64) 
yielding the form for the potential as 



ViA") - A^ [220 + 264C{-AY/^ + IGdC^i-A'f + C^-A'f/^ 



(65) 



where we have set an integration constant to zero since it would correspond to a 
cosmological constant. We have thus found that a time component of a vector field coupled 
to the Gauss-Bonnet invariant may generate 1) an early infiation epoch and 2) the ACDM 
background expansion without invocation of a cosmological constant or any scalar fields. 
This may happen when both the coupling and the potential are of a polynomial form, 
Eq. (63,65) in the case of dark energy and Eq.(60) in the case of inflation. These results 
seem quite interesting, keeping in mind that both vector flelds and the Gauss-Bonnet term 
have crucial roles in fundamental theories. Suppose a theory gives corrections to the action 
in a form of power-law expansion of the fleld strength \A\. Our results then imply that the 
negative powers could cause the early inflation and the positive powers could be responsible 
for the present acceleration of the universe. 
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3.3. Vortical perturbations 

The conventional decomposition of the metric separates perturbations into scalar, vector 
and tensor quantities according to their transformation properties under the spatial 
rotation group. Here we consider the vector-type (i.e. vortical) perturbations, and to 
avoid confusion we call them spin-1 perturbations. In fiat cosmology the metric, taking 
into account the spin-1 perturbations, can be specified without loss of generality by the 
line element 

ds"^ = a^{T) [-dr^ + B ■ rfxcir + cix ■ dx] . (66) 

We use here the conformal time t = J dt/a. Generally the spatial sections would also have 
inhomogcncities, but wc employ a gauge where they are fiat. The shift vector B then fully 
characterizes the inhomogeneities. It is transverse, V • --B = 0. The vector field we write 
as 

A^^a{T)[ip,L] (67) 

where the transverse vector L is a perturbation. The conformal background field </? = 
4>{T)/a{T) is introduced for convenience. It follows that 

and that = — v?^, which is understandable since is of spin-0 type. We can also write 
up to linear order that 

a'^R = 6{n' + a^R^B = 2A7fn', (69) 

where Ti = a' /a is the conformal Hubble factor and prime denotes derivative with respect 
to conformal time. The equation of motion for L then reads 



L" + 2nL' + 



e + {W + n^){l + 6a;') + ^f'n^n' + 2a V 



a2- 



0. (70) 



In this gauge the metric perturbations decouple from the evolution equation of the spin- 
1 type perturbation of the vector field. A wave-like propagation appears, unlike in the 
unit-norm models where the equation of motion instead becomes a first-order differential 
equation, whose solutions arc diluting [78]. Wc arc then able to deduce various properties 
of this perturbation. It is interesting to note that each component propagates like a scalar 
field fluctuation in the spatially flat gauge living in an effective potential given by the 
square bracket (excluding the gradient term). However, for the speciflc models considered 
in the last subsection the background equations motions reduce this to 

L" + 2nL' +[e + {n' + n^)]L = o. (71) 

Thus the evolution of the modes is independent of the detailed form of the coupling 
functions, and depends only on their effect to the background. As canonical scalar 
flelds, the components of L propagate with light speed at small scales. This important 
fact establishes that the spin-1 type perturbations are stable at small scales, quantum 
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mechanically consistent (i.e. not ghost modes^) and causal. For these models we find 
that the modes decay during inflation also outside the horizon, though at slower rate than 
vector modes usually, L ~ a~^'^^/^. However, under an expansion with We// > 1/3 
the superhorizon vector modes become tachyonic. This kind of effect could amplify the 
vorticity seeds during reheating. For radiation domination a constant mode exists. Given 
a regular primordial spectrum, vortical can result in large B-mode polarization signal at 
small scales £ ~ 2000, which is distinguishable from the effects of tensor and gravitational 
lensing on polarization [79] , and therefore possibly detectable with Planck. However, such 
effects might be pronounced in the present models. 

4. Conclusions and Discussion 

We investigated field theories where vector fields mediate non-gravitational interactions 
and drive the accelerated expansion of our universe, both during the early infiation and 
the present dark energy domination. The motivation for this investigation comes from the 
frequent appearance vector fields and possibile Lorentz violations in fundamental theories, 
from the need to the test the robustness of the basic assumptions of cosmology, and from 
the hints of statistical anisotropy of our universe that several observations seem to suggest. 

In this study we considered the question whether such vector fields could have an 
important cosmological role instead of the more popular hypothetical scalar fields and 
their exotic variants. A usual reason to exclude this otherwise reasonable possibility is 
that vector fields are incompatible with the established isotropy of the universe. However, 
space-like fields could feature only small anisotropy, which some anomalous observations 
indeed hint could exist in the universe. Secondly, time-like fields do not necessarily break 
isotropy, but they have non-trivial dynamics only if non-minimally interacting with gravity. 
To summarize our results, we found for space-like fields that 

• There exists for space-like fields (anisotropic) scaling solutions even without couplings 
to matter. 

• The cosmological bounds are satisfied in models of possibly large anisotropy if 1) the 
anisotropic era begins just recently or 2) the eccentricity cancels today. 

We learned about (time-hke) vector-tensor models that 

• With nonminimal couplings to gravity, dynamical solutions for cosmological vector 
fields become naturally available. 

• In particular, there are models with coupling to the Ricci or to the Gauss-Bonnet 
scalar that incorporate infiation and dark energy. 

In general, the amount of fine-tuning in vector dark energy models is similar to the 
amount of fine-tuning in scalar field models. As we have shown that there exists scaling 
attractors, the initial conditions for the fields and their velocities do not matter. One is 

% To see this, the equation of motion itself is not enough. One has to also check that the overall sign of 
the vector Lagrangian in the action corresponds to the right sign of the Hamiltonian, as now happens to 
be the case. 
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then left to explain the energy scales involved. Just like with scalar fields, by redefinitions 
and zero-point shiftings one may choose the numbers that go into the Lagrangian, but the 
field turns out to be extremely light. As an example, consider the exponential potential, 
V{A^) = Voe'^^^^^~^o\ If the potential energy is a significant contribution to the expansion 
rate today, the Priedmann equation tells that the potential is of the order V ~ H^/k. 
However, the potential scale Vq can be set to whatever seems natural by choosing suitably 
the parameter A^. The mass of the associated particle, ~ VV" ~ 10~^^ eV, if the slope 
A is roughly of order one. For the power-law potential, V{A'^) — Vo{kA'^)"', one finds the 
exact relation 

-— — 

ruA = (^YJ (72) 

where V is the value of the potential at a given time. If the time is nowadays and the scale 
of the potential is of order of the Planck mass, then the mass of the field is the Planck 
mass suppressed by a factor lO"^^^-^"-^/"'), which for large inverse powers can be quite tiny. 
These considerations hold also even if the potential vanishes [36], since one may interpret 
that the mass effectively arises from the nonminimal coupling. For the vector inflatons, 
one may find TeV scale masses and potentials. 

Thus, several cases of observationally viable and theoretically motivated vector fields 
in either FLRW or in a Bianchi I background were found. Those include models where 
again scaling solutions do exist for (space-like) vectors, with and without couplings to 
matter, and with given forms of the vector field potentials. With matter and the vector 
field present, there are three kinds of scaling solutions, which were summarized in Table (1). 
Allowing anisotropy of the dark energy equation of state could be useful in the quest for 
a realistic description of the present cosmological acceleration, since the exact symmetries 
of the FLRW metric exclude from the beginning for example these space-hke vector fields. 

We also considered vector fields non-minimally coupled to gravity. We showed that a 
time-like vector coupled to the Ricci scalar or to the Gauss-Bonnet invariant may generate 
both the infiation and the ACDM expansion. In the vector Gauss-Bonnet model, the 
form of the potential and the coupling which mimic the concordance model background 
expansion can be simple polynomial functions. Negative powers may drive infiation and 
positive powers the late time acceleration. Amusingly, both the time-like field and the 
Gauss-Bonnet term become trivial if the coupling is set to zero. Indeed, we have presented 
the first inflation models driven by time-hke vector fleld. This may have relevance to 
the problem of definition of the initial vacuum state of the universe. We considered the 
evolution of spin-1 type perturbations showing that it is both causal and stable. 
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